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Control of the transfer of quantum information encoded in quantum wavepackets moving along 
a spin chain is demonstrated. Specifically, based on a relationship with control in a paradigm of 
quantum chaos, it is shown that wavepackets with slow dispersion can automatically emerge from 
a class of initial superposition states involving only a few spins, and that arbitrary unspecified 
travelling wavepackets can be nondestructively stopped and later relaunched with perfection. The 
results establish an interesting application of quantum chaos studies in quantum information science. 
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I. INTRODUCTION 

Great effort is being devoted to studies of spin chains 
as promising "quantum wires" for quantum information 
transfer. With spin chains, a quantum state can be trans- 
ferred without requiring an interface between the commu- 
nication channel and a quantum computer i.e., quan- 
tum information can be transferred and processed with 
the same hardware. Spin chains also allow for quantum 
computing with an always-on interaction 0, 0] , even in 
the presence of a a global control field Q . The latest ex- 
perimental progress on fabrication and characterization 
of atomic spin chains was reported in Ref . Q . Spin chain 
Hamiltonians may be also realized by atomic gas in an 
optical lattice. 

Perfect state transfer in spin chains might occur un- 
der special circumstances [1,M[3- However, in the gen- 
eral case, dispersion effects often degrade the transmis- 
sion fidelity and improving the fidelity becomes a central 
issue. Notably, it has been proven that the transmis- 
sion fidelity can be significantly improved if the receiver 
stores the received signal in a large quantum memory 
before decoding [8,]. Another general approach to high- 
fidelity quantum state transfer advocates the use of quan- 
tum wavepackets to encode the quantum state of a qubit 
[9, 10]. This approach is important because dispersion of 
wavepackets can be insignificant. In particular, Osborne 
and Linden Q have shown that high transmission fidelity 
can be achieved by exploiting, if attainable, a Gaussian 
wavepacket whose shape is well preserved. The slow dis- 
persion of a wavepacket can be further suppressed by 
applying a static parabolic (hence global) magnetic field 

m 

In the context of the wavepacket approach to quantum 
state transfer, we focus below on two questions: (1) how 
can one create spin wavepackets with certain desired fea- 
tures, and (2) how can one control the motion of a quan- 
tum wavepacket in a spin chain so that the packet can 
be stopped at an arbitrary time, held, and then restarted 
later, without loss of quantum information. Such type of 
controlled quantum state transfer, if possible, should be 



a highly valuable tool in a variety of situations, e.g., cases 
in which the information receivers need additional wait- 
ing time to repair a quantum memory, or to prepare for 
a time window of high transmission fidelity. The impor- 
tance of stoppable quantum state transfer may be also 
appreciated by noting the analogy to the potential im- 
pact of the stopping of light [12] in quantum information 
science. Further, a working scenario for the stopping and 
perfect relaunching of quantum state evolution of a spin 
chain should be also of considerable interest in the con- 
text of perfect quantum state reconstruction and perfect 
quantum state storage in systems of interacting qubits 



In this paper we first show that by optimizing a 
particular transport property using quantum superpo- 
sition states comprising only a few spins (e.g., four or 
five), wavepacket pairs with some highly desired fea- 
tures emerge automatically from the ensuing dynam- 
ics. We then demonstrate that by applying a sequence 
of pulsed parabolic magnetic fields one can manipulate 
these wavepackets, stopping them and later relaunching 
the travelling wavepackets without individually address- 
ing the spins. As shown below, the stopping, followed by 
relaunching, can in principle perfectly preserve the quan- 
tum information being transferred. This is made possible 
by taking advantage of powerful relationships between 
controlling spin dynamics and controlling quantum dif- 
fusion dynamics in a paradigm of quantum chaos. 

This paper is organized as follows. In Sec. II we intro- 
duce a mapping between a Heisenberg spin chain kicked 
by a parabolic magnetic field and a paradigm in quan- 
tum chaos In Sec. Ill we propose a conceptually 
simple approach to the creation of spin wavepacket pairs 
moving along the spin chain with slow dispersion and 
other desired features. The key result of this work is in 
Sec. IV, where stopping and relaunching spin wavepack- 
ets are studied both numerically and analytically. Section 
V concludes this paper. 
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II. HEISENBERG SPIN CHAIN IN A PULSED 
MAGNETIC FIELD AND THE DETLA-KICKED 
ROTOR 

Consider then an open-ended Heisenberg chain of N 
spins in a constant magnetic field B and subject to a 
paraboUc (5-pulsed magnetic field. The Hamiltonian is 
given by 

J N-l N 
n—1 n— 1 

3 n=l 

where a = (u^, cr^, fi^) are the Pauli matrices, J is 
the nearest- neighbor spin-spin interaction constant, Cj 
and riQ are the coefficient and minimum location of the 
parabolic kicking field, and Tq is the kicking period. Be- 
low, we denote the n = 1 {n = N) spin as the left (right) 
end of the chain. The constant field B lifts the system 
degeneracy and the dynamics is restricted to a subspace 
with fixed total polarization Sz defined as 

N 
n=l 

Throughout this work we consider only the subspace of 
= 1. 

Let |m) be one of the basis states, with the mth spin 
up and all other spins down. The propagator for the time 
period [jTo - 0+, (j + l)To - 0+] is V{2JTo)U{Cj). Here 
U{Cj) represents the action due to the delta pulse, with 

(m|;7(C,)|n) = exphi(Cj/2)(n - nof]S,nn. (3) 

The term V stems from the evolution inherent in the 
Heisenberg interaction. An important recent study [f3| 
has shown that, in the N +oo limit (and apart from 
some irrelevant phase) 

(m|l/(2JTo)|n) « J(„_„) (2 JTo), (4) 

where J(m-n) is an ordinary Bessel function. The ana- 
lytical behavior of U{Cj) and V{2JTo) is therefore com- 
pletely in parallel with that associated with the propaga- 
tor of the (5-kicked rotor (DKR) (the best known model 
in quantum chaos [lB|) with Hamiltonian 

Hdkr = {P- P^fl2 - K cos{e) J2 ^(t - j)- (5) 

i 

Indeed, in the representation of the basis states |m) = 
cos{m6)/y/Tr and for an effective Planck constant h, the 
DKR propagator takes the familiar form v{k)u{fi), with 

{m\u{h)\n) = exp[-i(n/2)(n - haf]5mn, (6) 



and 

(m|{;(fc)|n) «z(™-")j(,„_„)(fc) (7) 

with {k = K/h). Comparing these two systems, it is 
clear that upon the mapping 

|m) ^ |to), (8) 

2JTo ^ fc, (9) 

no ^ Po/h, (10) 

^ h, (11) 

the many-bod y sp in chain dynamics is mapped to that of 
DKR p^. [Tsl. fl7|. i.e., the motion of a spin wavepacket 
along the spin chain is mapped to DKR quantum dif- 
fusion dynamics in its m-space. Hence we can, when- 
ever possible, shed light on the former by considering as- 
pects of the latter, e.g., quantum resonance, Kolmogorov- 
Arnold-Moser (KAM) curves in phase space, etc. More 
significantly for this work, as shown below, it allows us 
to use tools from the control of quantum DKR dynamics 
[Tsl [T9I [20j to manipulate states in the spin chain. Fur- 
ther, this mapping between spin-chain and DKR allows 
us to go beyond the parameter regime confined by the 
true DKR (discussed below). 

III. GENERATION OF SPIN WAVEPACKETS 

In the context of the wavepacket approach to quantum 
state transfer, we now consider the first issue on spin 
wavepacket generation [21|. Given a small number of 
basis states that could be used for encoding the state 
of a qubit, what initial superposition states should be 
exploited to induce the creation of quantum wavepackets 
with slow dispersion? Here this interesting question is 
considered in the absence of an external field, where the 
system propagator is given by V{2JTo). Remarkably, the 
associated DKR analogy now becomes a case of quantum 
resonance with h = ^tt, with a propagator analogously 
given by 

v{k = 2 JTo) = exp[i(2jro) cos(6i)]. (12) 

Using this connection, the issue becomes to find initial 
superposition states within a given small subspace, such 
that the evolving quantum state remains well localized. 
At first glance this "localization" requirement seems too 
demanding because the main feature of quantum res- 
onance dynamics is ballistic diffusion in the DKR m- 
space. However, as we have discovered, this can still be 
obtained by maximizing a diffusion rate of DKR. Qual- 
itatively speaking, for superposition states maximizing 
a quadratic diffusion rate introduced below, the ensu- 
ing dynamics will push outwards as much as possible the 
excitation profile in the m-space, thus generating two 
well-separated wavepackets with almost zero amplitude 
in between. 
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FIG. 1: Emergence of a wavepacket pair in a Heisenberg chain 
of 201 spins, shown with the projection probabihty of the 
many-body wavefunction onto basis states |m). The initial 
condition is a superposition state \^mo) given by Eq. (15), for 
mo = 101 in (a) and mo = 30 in (b). Th system wavefunction 
then evolves, with its shape given by the solid lines at time ti 
with 2Jti = 15, and by the dashed lines after an additional 
period t2, with 2Jt2 = 30 in (a) and 2Jt2 = 45 in (b). The 
arrows show the travel direction of the wavepackets. 



Quantitatively, let us first define the diffusion rate op- 
erator as 



(13) 



where E{t) is the energy operator for the free rotor in the 
Heisenberg representation. For the quantum resonance 
case considered here one obtains 



D = Asm^{9), 



(14) 



where A is a constant. Note that D only couples states 
\m) of the same parity. Consider now a sample case where 
a superposition state 



n=2 



P2n\mo + 2n) 



(15) 



is exploited to encode a qubit state, i.e., only five basis 
states are used here. The state with the largest diffusion 
rate, denoted D, is simply given by the eigenfunction of 
D in the subspace of |mo -|- 2n) (—2 < n < 2) with the 
largest eigenvalue. In particular, if these basis states do 
not involve state |0), then the maximized D is attained 
if Po = 0.577, /3_2 = P+2 = -0.5, /3_4 = P+a = 0.289. 

The significance of such an initial superposition state 
with maximized D is demonstrated in Fig. 1(a), with 
niQ = 101 for a 201-spin chain. In particular, a well- 
separated wavepacket pair is seen to quickly emerge, and 
its dispersion after its emergence is impressively slow in 
the absence of any external static fields. Note that, unlike 
the accelerator mode approach proposed in Ref. [ij , the 
wavepacket pair is created by the system dynamics itself. 
Note also that the excitation amplitudes between the two 
wavepackets are surprisingly small. Because the total po- 
larization here is fixed, a well-separated wavepacket pair 



directly indicate quantum entanglement between well- 
separated parts of the spin chain, and their travel in 
opposite directions distributes information or entangle- 
ment to both ends of the spin chain. Further results 
(not shown) indicate that in the case of minimized D or 
an arbitrarily chosen initial state, the system dynamics 
generically creates a quickly delocalizing state (note also 
that even in chaotic cases different initial superposition 
states may also lead to dramatic differences in the en- 
suing quantum diffusion dynamics (2^). These further 
demonstrate the important role of an initial superposi- 
tion state with a maximized diffusion rate in encoding 
the quantum information. Certainly, if more basis states 
are allowed in encoding the quantum information, then 
wavepacket pairs with even slower dispersion can be cre- 
ated with the same approach. 

It is also desirable to be able to create a well-separated 
wavepacket pair that transfers information to a common 
end of the spin chain. For example, if two wavepackets 
with identical shape can be created, then one of them 
may be analogous to a "backup" copy as the other is be- 
ing transferred and received first. Note that this possibil- 
ity is not in violation of the quantum no-cloning theorem, 
because here the two wavepackets do not independently 
describe the quantum state of the involved spins. Rather, 
the two wavepackets describe the entanglement between 
two particular sections of the spin chain. 

The creation of such a wavepacket pair is achieved here 
by going beyond the kicked rotor perspective and exploit- 
ing the boundary effect associated with the spin chain. 
That is, we apply the above scenario, but with the ini- 
tial encoding state \ipmo) located at mp < N/2, and with 
the requirement that no information receiver presents at 
the left (n = 1) end. To be more specific, consider a 
sample result shown in Fig. 1(b), with mg = 30. The 
wavepacket creation dynamics in the early stage is seen to 
be analogous to the case of Fig. 1(a). Sometime later, the 
generated wavepacket moving to the left hits the bound- 
ary and gets reflected. As demonstrated in Fig. 1(b), 
this then creates a pair of wavepackets where both of the 
members of the pair are moving to the right, with their 
shape indistinguishable from one another, with almost 
zero excitation in-between, and a peak-to-peak distance 
of 2(r7io — 1) spins. 



IV. STOPPING AND RELAUNCHING 
WAVEPACKET PROPAGATION 

In this section, taking the wavepackets obtained in the 
previous section as examples, we shall consider an in- 
dependent issue, namely, stopping and relaunching the 
quantum state transfer along a spin chain. A parabolic 
magnetic field [see Eq. ([T|)] is proposed as the control 
field with a simple global feature, and we aim to achieve 
our control objective with an always-on Heisenberg in- 
teraction. 

Thanks to the DKR analogy discussed above, the prob- 



4 



-g 0.08 
Q. 0.04 



-g 0.08 
0.04 



.g 0.08 
Q. 0.04 



- 


k k 


(a) _ 




1 


51 101 151 


201 






(b) _ 




1 


51 101 151 


201 




1 1 1 

, ^. .JVimA^. .-^ , 


(c) _ 





1 51 101 151 201 

index of spins 



FIG. 2: Stopping of quantum state transfer in a spin chain by 
a parabolic kicking field, with Cj = 0.5. The initial state, the 
meaning of the wavefunction profile, and the spin chain used 
for calculations are the same as in Fig. 1(a). 2JTo — 0.25, 
and the first 5-kick comes at time ti with 2Jti = 15. (a) and 
(b) are the results after 100 and 200 kicks. The kicking field is 
then turned off and (c) displays the state after an additional 
period t2 with 2Jt2 ~ 30. 



lem is now converted into the question of how to stop and 
successfully restart the transport process in the DKR m- 
space. We do so below by exploiting control features of 
the DKR system. Two features are relevant, one classical 
and one quantal. 

The classical basis of our control scenario arises by 
exploiting the phase space KAM curves of the underly- 
ing classical dynamics. That is, if the kicking magnetic 
field is sufficiently frequent that the chaoticity parame- 
ter |2 JToCj I in the classical DKR is sufficiently small, the 
underlying classical dynamics will be mainly integrable 
and the associated KAM curves will present strong bar- 
riers to the quantum transport in the m-space. Because 
KAM curves will be almost everywhere, these classical 
structures can effectively stop the travel of arbitrary and 
unknown quantum wavepackets. 

Figure 2 displays the fate of the moving wavepacket 
pair shown in Fig. 1(a) (solid lines) after a kicking 
parabolic field is introduced. As is clearly seen in Fig. 
2(a) and Fig. 2(b), the transfer of the wavepacket pair 
to both ends of the spin chain is stopped. This dynam- 
ical effect can also be understood as a type of quantum 
Zeno effect achieved by frequently applying (but far from 
infinitely fast) external pulses to a system. 

However, although the wavepacket pair in Figs 2(a) 
and 2(b) has stopped moving, its internal structure is 
seen to be changing in a subtle manner. This indicates 
that evolution of the quantum phases characterizing the 
stopped wavepackets is still not frozen. This fact turns 
out to be disastrous when the kicking field is turned off 
in order to relaunch the state transfer. For example. Fig. 
2(c) displays the wavefunction after the kicking field has 
been off for a period of t2 with 2 Jt2 = 30: the background 
fluctuation is greatly increased, and the main peaks of the 
wavepacket pair do not move further. 

Hence, using KAM curves alone, which is a purely 



TABLE I: The j-dependence of Cj [see Eq. (1)] in an ex- 
plicitly designed pulse sequence for the stopping of arbitrary 
wavepackets for a period of 2M kicks. C is a constant dis- 
cussed in the text. Note that some system parameters used 
here are beyond what is allowed in a true kicked rotor system. 
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classical control mechanism, does not offer a satisfactory 
means of stopping the wavepacket. To improve the con- 
trol one must compensate for the quantum phases that 
are accumulated during the stopping process. This phase 
accumulation is due to the spin-spin interaction as well 
as the kicking field. 

Consider then an important observation made in our 
previous work on the quantum control of DKR dynamics 
til,|2Q|, I.e., 

1 /-^^ 

{■m\v{k)\n) — — cos{m9) exp[ik cos{9)] cos{n6)d9 

TT Jo 

= (m|(-l)™z)(-fc)(-l)"|n) 

= {m\u{2TT)v{~k)u{2TT)\n). (16) 

The first line in Eq. holds by definition, the sec- 

ond line becomes obvious if the integration variable 9 is 
changed to 9 + tt, and the last line is obtained by use of 
the definition of u(h). Equation (jl6p hence proves 

v{k) = u(27r)i)(-fc)u(27r). (17) 

Returning to a finite spin chain system, this result indi- 
cates that 

V(2JTo) w U{2Tr)V{-2JTo)U{2Tr). (18) 

That is, the sign of the intrinsic interaction constant J 
can be effectively reversed if we apply two parabolic i5- 
kicks of particular strength. As such, it becomes possible 
to compensate for the quantum phase evolution inherent 
in the spin chain. As to the quantum phases induced by 
the kicking field, they can also be compensated for by 
considering kicking fields with the sign of Cj reversed. 

Given these considerations we present in Table I an 
explicitly designed special pulse sequence that can re- 
launch stopped wavepackets with perfection. For this 
special pulse sequence, the KAM curves associated with 
small \2JTqC\ still play a key role because they directly 
prevent the state transfer, in the same manner as demon- 
strated in Fig. 2. What is remarkable now is the total 
time evolution operator associated with the entire stop- 
ping process. In terms of the DKR analogy, this operator 
can be written as (after some manipulation) 

[^-C/2)v{k)ui-C/2)f' [uiC'/2)vik)u{C'/2)f^ 
= [ui-C/2)vik)u{-C/2)]^' [u{C/2)vi-k)u{C/2)]^' 
= [ui-C/2)v{k)u{-C/2)]^'' [u{C/2)v{-k)u{C/2)f'' 
= ••• = 1, (19) 
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FIG. 3: Same as in Fig. 2, but for an explicitly designed 
pulse sequence given in Table I, with C = 0.5, M — 100. The 
transfer of the wavepacket pair along the spin chain is stopped 
for 2A1 kicks. Panels (a) and (b) are the results after 100 and 
200 kicks. The result in (c) shows that the state transfer is 
indeed relaunched with perfection. 

where C' = C + An, M' = M - 1. In obtaining Eq. HH) 
we have used 

u{C) ^u{C/2 + 2tt)u{C/2 + 2tt) (20) 

and Eq. ^TE\\ . Equation (fTO|) proves that at the end of the 
stopping time all properties characteristic of an unknown 
quantum wavepacket can be exactly restored. This exact 
rephasing indicates that the dynamical evolution associ- 
ated with the second AI/2 kicks, in addition to offering 
a dynamical barrier to stop the quantum transport, pre- 
cisely reverses the evolution associated with the first M/2 
kicks. As such, the stopping is entirely nondestructive, 
as long as the system is not subject to noise effects dur- 
ing the stopping process. Evidently then, wavepacket- 
assisted information transfer can be perfectly relaunched 
as the kicking field is turned off. This theoretical result 
applies exactly to an infinitely long spin chain. But fortu- 
nately, as also demonstrated below, it applies extremely 
well to a finite-length chain. Note also that the designed 
pulse sequence in Table I is a significant extension beyond 
a true DKR system because both positive and negative 
"h" {h ^ Cj) are exploited here. 

In parallel with Fig. 2, Fig. 3 displays a computational 
example using the pulse sequence given in Table I. As in 
Fig. 2, right before the first kick the quantum state is the 
wavepacket pair described by the sofid lines in Fig. 1(a). 
Figure 3(a) confirms that the wavepacket pair is not mov- 
ing as the special pulse sequence is on. Figure 3(b) shows 
that at the end of the stopping period, we have restored 
the initial condition [compare Fig. 3(b) with the solid 
lines in Fig. 1(a)]. The restoration fidelity in the numer- 
ical calculations for a 201-spin chain is found to be higher 
than [1 — 10"^"^]. The kicking field is then turned off. As 
expected, quantum state transfer is relaunched and the 
wavepackets continue their journey, with slow dispersion, 
towards both ends of the spin chain [Fig. 3(c)]. Indeed, 
results in Fig. 3(c) are indistinguishable from the dashed 
lines in Fig. 1(a). 



The control scenario proposed in this work can also 
lead to other very interesting approaches to the manipu- 
lation of quantum entanglement dynamics of a spin chain. 
Here we briefly discuss three possibilities. First, by mod- 
ifying the kicking field profile we can choose to stop only 
one component of a wavepacket pair, e.g., of the pair 
shown in Fig. 1(b) with dashed lines, thereby offering an 
interesting method of tuning the time delay between two 
wavepackets moving in the same direction. This then 
offers a means of controlling the distance between two 
entangled parts of the spin chain. Second, because the 
sign of the intrinsic spin-spin interaction constant J can 
be effectively reversed if we apply (5-kicks of particular 
strength, it can be easily shown that one can bounce 
back an arbitrary and unknown moving wavepacket to 
the sender at a time of our choosing. Third, by control- 
ling the time delay between the two wavepackets and/or 
taking advantage of the feasibility of time-reversal, we 
may also recombine two localized wavepackets at a lo- 
cation different from that of the initial state. This re- 
combination dynamics resembles that of a double-slit ex- 
periment, thereby generating interesting interference pat- 
terns along the spin chain. Such kind of interference pat- 
terns of spin excitations, and their fate under a variety 
of circumstances, may work as a novel interferometer for 
fundamental studies in quantum physics. 



V. CONCLUSION 

To conclude, based on a mapping between a kicked 
spin chain and the delta-kicked rotor system [13], we 
have shown that previous quantum control results in the 
delta-kicked rotor system [l^ El) HO, HI] can be applied 
to the control of spin wavepacket propagation and hence 
the control of propagating quantum information encoded 
in wavepackets. Specifically, we have proposed a simple 
approach to wavepacket creation in a Heisenberg spin 
chain and demonstrated the possibility of stopping and 
relaunching information transfer without individually ad- 
dressing spins or turning off spin-spin interactions. Sev- 
eral interesting applications of this work in manipulating 
the dynamics of a spin chain are also discussed. The re- 
sults indicate that many insights from the quantum chaos 
research can be very useful for quantum information sci- 
ence. This work also adds more support to the use of 
spin chains as quantum wires, and might be useful in 
designing new quantum computation algorithms with an 
always-on qubit-qubit interaction 0, y]. Extensions to 
other types of spin chains are under consideration. 
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